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Perturbation theory of the dynamic inverse spin Hall effect with charge 

conservation 
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We present gauge-invariant theory of the dynamic inverse spin Hall effect driven by the 
spin-orbit interaction in metallic systems. Charge conservation is imposed diagrammatically 
by including vertex corrections. We show the charge current is induced by an effective electric 
field that is proportional to the spin current pumped by the magnetization dynamics. The 
result is consistent with recent experiments. 
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1. Introduction 

Control and manipulation of spin current, the flow of the spin degree of freedom in solids, 
has been extensively studied since the discovery of the spin Hall effect in semiconductor 
systems. The spin Hall effect converts electric voltage into spin current via the relativistic 
spin-orbit interaction. The induced spin current was first detected indirectly by observing 
spin accumulation at the edge of the sample.^' ^ The spin-orbit interaction also works in the 
opposite way; that is, when the spin-orbit interaction acts on the spin current, charge current 
is expect to arise. This is the inverse spin Hall effect. By use of the inverse spin Hall effect, 
direct electrical detection of the spin current has become possible.^' ^"^^ 

Depending on the creation mechanism of the spin current, there are two inverse spin 
Hall effects. One is the static inverse spin Hall effect involving non-local spin injection. ^"^'^ 
Spin current is created from the electric field via the direct spin Hall effect and current 
arising from the inverse spin Hall effect can be observed by attaching a nonmagnetic contact 
of Pt or Au having strong spin-orbit interaction. The other effect is the dynamic inverse 
spin Hall effect involving a spin pumping mechanism^^"^^ with dynamic magnetization. The 
concept was employed by Saitoh who used a microwave to excite the uniform precession of 
the magnetization. The resulting electric voltage was found to be perpendicular to the spin 
current.^' 

The creation of electric current by dynamic magnetization is well known in classical elec- 
tromagnetism as Faraday's induction law. In contrast to Faraday's law, which is classical and 
macroscopic, the inverse spin Hall effect is a quantum relativistic effect that is induced in 
solids. The crucial practical difference between the two is that while the current induced by 
Faraday's law diminishes when the system size is reduced with the field kept constant, the 
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efficiency of tlie inverse spin Hall mechanism remains constant. The use of the dynamic inverse 
spin Hall effect is therefore expected to be essential in high-density integrated devices. 

Theoretical study of non-local spin injection was carried out by Takahashi and Maekawa.^'^ 
Their discussion was based on a semiclassical argument that the charge current and spin 
current are related by the spin-orbit interaction j oc cr x j's, where cr and jig represent the 
spin polarization and flow direction of the spin current, respectively. Their work qualitatively 
explained the experimentally observed static inverse spin Hall effect. However, the theory, 
which simply predicts the generation of charge current when spin current and a spin-orbit 
interaction exist, violates the conservation law for charge. Charge conservation is a crucial 
issue in any physical phenomenon and a framework that upholds the conservation law is 
certainly required. 

The dynamic inverse spin Hall effect was studied theoretically for two-dimensional dis- 
ordered electron systems with Rashba interaction.^^' Coupling to the localized spins was 
treated perturbatively assuming weak coupling as is the case for ferromagnet/normal metal 
junctions. The dominant contribution to the current was shown to be Jq, ~ \^^a(izk[.S x S')^)^ 
where S denotes a local spin. Ar and ea^z are the strength of the Rashba interaction and 
anti-symmetric tensor respectively, ( ) represents the average owing to electron's diffusive mo- 
tion and the z axis is perpendicular to the plane of the junction. The result seems reasonable 
since the current arises from the energy dissipation (damping) of the spin system, S x S. 
The obtained current direction is consistent with experimental results.^' It is, however, an 
open and critical issue whether the Rashba inverse spin Hall mechanism applies to the ex- 
perimental systems,^' since the Rashba interaction known so far in metallic systems occurs 
on the surface of nonmagnetic metals^^ and not at the interface between a ferromagnet and 
nonmagnet. 

Gauge-invariant (charge-conserved) theory for the dynamic inverse spin Hall effect in 
metallic (i.e., symmetric) spin-orbit systems is therefore necessary. The aim of this paper is 
to present such theory on the basis of microscopic quantum many-body theory. Charge con- 
servation is guaranteed by evaluating Feynman diagrams while keeping the gauge invariance. 
This means we include vertex corrections, which represent the effect of electron diffusion. As 
the spin-orbit interaction, we consider that induced by random impurities, which is dominant 
in metallic systems. A crucial feature of the present system is that the translational and rota- 
tional invariances are recovered after averaging over the random impurities, in sharp contrast 
to the case for a Rashba system with explicitly broken symmetry. ""^^'"^^ 

We will show that the induced current is written as 

ja = aisneais-yj^^fs - DdaP, (1) 

where jj^ is the spin current density locally created by the dynamic magnetization, aisH 
represents the efficiency of the spin-charge conversion, D = is the diffusion constant 
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and p is the charge density, (ep and r the Fermi energy and the hfetime, respectively) The 
coefficient aisH and expressions for jj^ and p are given below, eq. (1) is exact within second 
order perturbation theory and for slowly varying spin structures. The result satisfies the charge 
conservation law p + V • jf = 0. The first term of eq. (1) clearly indicates that spin-charge 
conversion indeed occurs in dynamic inverse spin Hall systems. This spin-pumped part is 
local and thus the spin current acts as an effective electric field. In contrast, the spin current 
in static inverse spin Hall systems arises diffusively from spin accumulation^''^ and thus the 
static inverse spin Hall current is expressed by a gradient term Vs, where s is the electron 
spin density. In the present dynamic case, diffusive spin current, proportional to Vs, does not 
exist. The physical origins of the spin currents for the dynamic and static inverse spin Hall 
effects are therefore different. 

2. Model 

The system we consider is a disordered electron system with lifetime r in the presence 
of a random spin-orbit interaction and coupled to dynamic localized spin S{x,t). The total 
Hamiltonian is given by 

where tp^ and ■0 are the electron creation and annihilation operators. The first term describes 
conduction electron and -ffimp = / d^x'ip{x,t)^U{x)Tp{x,t) represents elastic impurity scatter- 
ing, where U{x) is the random impurity potential. We assume averaging over the impurity 
position is carried out in the standard manner as {U{x)U {x'))ij^p = UQUimp^ix — x'), where 
uq and nimp are the strength of the impurity potential and the impurity concentration, respec- 
tively. This impurity scattering f^imp gives lifetime as r = 2-KN{Q)n- — i?' '^'^^'^^ -^(0) = i^J^ 
is electron's density of states at the Fermi energy and V is the system volume. The electron's 
diffusion motion expressed by the diffusion ladder also arises from the impurity scattering. 
The coupling to the localized spin is represented by the s-d type exchange interaction. 



^^cx = Jcy. j d^xS{x,t) ■ '4)\x,t)(T^{x,t). 



We assume this interaction is weak (compared with the inverse lifetime), ^^^^ ^ 1. In the 
case of ferromagnet-nonmagnet junctions, the electrons contributing to the measured current 
are those in the nonmagnet, and the exchange interaction acts only at the interface. The 
assumption of a weak exchange interaction would thus be correct. The spin-orbit interaction 
we consider is due to random impurities: ^^"^'^ 

^so = -^Ciifc j d^xdiU{x) (tl;\x,t)a''djip{x,t)'^ , 

where ip"^ djil^ = ip'^djip — ((9j'(/'^)^, Ago represents the spin-orbit strength and U{x) is same 
impurity potential as in //imp- The spin-orbit interaction is treated to the first order. 
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3. Method 

Let us calculate the electric current induced from the magnetization dynamics. We assume 
a slowly varying spin structure in both space and time, and expand the current with respect 
to dt and V. The total electric current is given as ja = ja + Ja°) where 

2m 



J. 



SO 



eAs 



'■eaim.{ipH'^^d„iU)ip), 



h 

and the charge density as p = e(^^'0). Here, ja is the normal term (skew-scattering term) and 
Jq*-* is the anomalous term which arises from the spin-orbit interaction due to the impurity 
(side jump term). The expectation values are written by the lesser Green's function as 

k,q 

fj'ix,t) = ^Y1 ^^^^^rna {k' - A:)' lim(yfc,_fctr[<T^G<_^ ,^,^^(t,t')])impe-^^''" 

kk'q 

p{x,t) = -i^^(trG<_^ ,,^^(t,0)impe-^''-". 

k,q 

Here, trace (tr) is over spin, and (. . .)imp represents the impurity averaging. Lesser Green's 
function G^y,i{t,t') is calculated by use of the contour ordered Green's function Gk,^y{t,t') = 
jfi{Tc[ck{t)(^]^i{t')\) , where Tc is the contour-ordering operator along the contour C of the 
complex time.^^ The lesser Green's function satisfies the Dyson given as 

Glu'M =&-t')h,k' 

+ Jo.Y.\l "^^^Skit - ti)5;_fc(ti)a^Gp,fc.(ii,t') 

+ ^ f^fc-p(^ P)i (yjdti9k{t - h)a'Gp,k'ih,t')j , 

where gk{t — t') is the free Green's function including the lifetime. The Fourier component 
of the impurity potential is defined as U^'-k = / dxU {x)e~'^^^~^'^'^ . In the Fourier space 
(defined as gk{t — t') = J ^gk,u}^~^'^^^~^ ^)i the free lesser component is obtained as g^^ = 
fi^)i9k,uj - 9k,J^ where f{uj) is the Fermi distribution function, and gl^^ and g%^^ = {gl J* 
are the retarded and the advance Green's function respectively. The retarded Green's function 
is given as g]^^ = [i] = ^). (The free Green's functions are diagonal in spin 

space.) 

3.1 Linear order in J^x 

The contributions to the charge current at the linear order in both Jcx and Ago are shown 
in Fig. 1. These contributions arise from the second order in impurity potential U^-k' (One 
of them is arising from Hgo, the other one from i^imp)- The vertex corrections to the normal 
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q,!2 

k+ql2,oH-Q 



Fig. 1. Diagrammatic representation of the charge current at linear order in the exchange interaction. 
Crosses and filled circles represent the current [j) and spin-orbit interaction (SO), respectively. 
The intersection with a filled circle represents the contribution of the anomalous velocity (j^'-') 
from the spin-orbit interaction. Solid and wavy lines denote the impurity-averaged electron Green's 
function and the exchange interaction S, respectively. Dashed lines denote the impurity averaging 
between the spin-orbit interaction and the impurity scattering. Complex conjugate to each process 
(a) and (b) are also included in eq. (2). 



charge current at the linear order in Jcx vanish exactly. The charge current at the hnear order 
in the exchange coupling Jcx is thus obtained as 



(2) 



where Sq{Q) is Fourier component of the localized spin, defined as S{x, t) = 
Q S'q(0)e*(^*~'^'^\ Contribution from an each diagram (a) and (b) of Fig. 1 are given 

as 



C 



1(a) 

ajl3;q,n 



m 



{k' - ky / (^j ^gk^l^u;9k'~i,u;9k' + ^,u;+n9k+i, 



C 



m 

ajl3;q,n 



k,k' 
duo 



-^9k~^,oj9k'~^,ui9k'+^,uj+n ) + ( / ■:^9k~^,ui9k+ 



duj 



2ir- 



^'Q9k'+%,Lu+n 



aji3-qn contribution of the normal charge current (ja), and ^^jl-qu is arising from the 

additional charge current (Jq*^). Using Langreth theorem, we obtain 

< r rh,, r 



'2^9k~ f ,uj9k'~^ ,uj9k'+^ ,uj+n9k+ § 



duj 
2^ L 



+ fi^)9l-i^^9l'^!i^^9l'+i^^+n9l+i^^+ 

(3) 

Since the contribution containing only or only is small (higher order of after inte- 
grating over k and k', the dominant contribution of Clj'j^g.q q come from the second term of the 
right-hand side of eq. (3). By using the approximation /(a; -|- $7) — f{uj) ~ Q,f'{io) (justified if 
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Q <^ u) and f'{uj) ~ —5{uj) at low temperature, we obtain 

k' 

where 5^ = gl. etc. The other contribution C^j^.^ ^ is similarly obtained as 



k,k' 



2^^i°'^^^''^'^^' ~ ^) ■ ^9k-i9k'-i9k'+i9l+i 

k,k' 



Here, we used the identity (7^_a — 9I.+ Q = (^^fc-| ~ ^fc+f ~ ^^'^j Sk-S-^k+S-' assume the 
smooth magnetization profile {\q\l <C 1, where / is the electron mean free path), and expand 
^aj^3-q n ^^^^ ^ai/3 <j H ^i^h respcct to wave vector q. We also assume that the system is spatially 
symmetric. We then obtain 

Ci';^.,,,n=Cl%{q,n) = ^S^,q^Y.'k\9i\W. (4) 

fcfc' 

Therefore the contribution of the skew-scattering and side jump are same in the slowly varying 
limit and at the linear order in J^x- Integration over k, k' is carried out using ^fcfifclfl'fcP = 
-^(0)^^! Y^k' bfc'P = iyMZL. "vVe finally obtain the current eq. (2) as 

j^^\q, n) = aisnce^p-ynq^Slin), (5) 

where c = g^V and oish = 

3.2 Second order in Jex 

We next calculate the contribution at the second order in the exchange coupling Jex. Total 
charge current at the second order in Jex is given as 

(6) 

(2) 

Here the contribution of P„/„no o is the local contribution which does not contain the 

vertex corrections as diagrammatically shown in Fig. 2. The diffusive contribution of the 

(3) 

charge current is presented by ^^^/g q r2i ' '^'^i^^ contains the vertex corrections to ja , and 
are shown in Fig. 3. 
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k'+q/2,(o+Qi 
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k'+Q/2,(o+Q 

Fig. 2. Diagrammatic representation of the charge current at the second order in the exchange inter- 
action. Complex conjugate processes are also included in eq. (6). 

3.2.1 Local contribution 

We calculate first the local contribution of the charge current. This contribution is given 
as 7^(2) =-d2(«) ,^2(fe+c) .^2(d) , 



^ai;q,Q,ni,n2 ~ ^ 

k,k' 



2(b+c) 

ai\q,Q,Q,\,Q.2 



m 



k,k' 

dio 
2^ 



9k^a.^uj9k+^,io+ni9k'+^,w+Qi9k'+a.^uj+n9]^+s.,uj+n 



k,k' 



duj 



a ,<^5fc/_ a ,(^5fc' + f ,w+Qi 5fc'+f 5k+Q ,a;+n 

'^ai!l,Q,ni,n2 = ^ihn^OLm ^[k' -k) (j — C/fc^fi^^^fc+f ,a;+ni5fc+a,^+nf fc' + f 

+ (/ -^9k^^^^9k'^^,u,9k'+l,u+m9k'+^,^+(!^ > 

where Q. = VLi + Q.2. Each contribution corresponds to the contribution labelled in Fig. 2. 
After taking the lesser component, we obtain the leading contributions as 



V 



2(a) 

ai;q,Q,V.\,Vl2 



2'Km 



fe.fc' 
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Fig. 3. Diagrammatic representation of the charge current with vertex correction due to the electron 
diffusion ladder (C) shown by double dotted lines (shown in (a),(b),(c)). There are two diagrams of 
each type (a), (b) and (c) due to the complex conjugate process. The diffusion ladder ( ( shown in 
(d) ) describes the diffusive electron motion and results in non-local contributions, which sum to be 
DdaP- These processes containing ( are the vertex corrections that are essential for guaranteeing 
the charge conservation. 



27r 

k,k' 

where c.c denotes the complex conjugate. After expanding with respect to the wave vector q 
and Q (assuming \q\l, \Q\l <C 1), we obtain 

fc,fc' 



^Qi;q,Q,ni,n2-*^^^i^«"™ p Z^^k\gk\ \9k'\ + ( 2 Z^^i^^dkl \9k' 

\ k,k' k,k' 
k,k' 

Using ^fc fc, efcbLPbL'l^ = Ek,k' efcbfcl^l^L'P = ^^^^^^^ we obtain 



,2{a) , ^2{6+c) _^2{d) _ .iV(0) TTep , 



We see that skew—scattering contribution (^Ji-^q^Q^fi-^^Q, 

ai;q^Q,r2i,Q2) ^'^^ ^^'^^ jump 

contribution (^^^^^ q ^^^^^ J are same. 
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3.2.2 Diffusive contribution 

We evaluate here the contribution which include the vertex corrections as shown in Fig. 
3(a)(b)(c). This contribution T^ai gQ^i ^2 given as 



V 



(3) 



1^=1 kk' 



(8) 



^lm,3{a) . ^lm,3{b+c) 



where 7 = I^a^^ /.^ = ^ ^nd 



-,Zm,3(fe+c) _ 
^q,Q,LO, 0.1,0,2 / J 



><9k,^3^,u9k,+^,Lu+nj9k',+^,^+f^i9ki+S.,u+n9k,+q.,u+n 



kv kj^ 



X9k,-.S.^^9k'^^S.^^9ki+^,uj+ni9k,+^,iu+ni9k,+q.^^+n- 

An infinite series over the index v gives the vertex corrections shown diagrammatically in Fig. 

(3) 

3(d). After calculating the lesser component of the T^ai-qQUi 02' obtain 



V 



(3) 

ai;q,Q,Ui,U2 



m 



u=l kk' 
diO 

2^ 



-,lm,3{a),< 



lm,3{b+c),< 



n: 



n=l9k^_„+^,uj+n 



Here 

-,lm,3{a),< 



kjy k' 



-,im,3(b+c),< 
^q,Q,Oi,n2 



kf ^k'l/ 



[-'2k'ik:; + [k',-k,)' 



Q-q 



9k,_3.9ki-^9ki+^9k,+!i9k^^ 



Expanding Q'^^qq^'q^ and Q'^qq^q^^ with respect to Q and q, we obtain 



-,lm,3{a),< 
^q,Q,ni,n2 



-Q"^q^6iu;)n,Yek'j9k\ 



ku kj^ 



9k' 



9k', 



and 



plm,3{b+c),< 
'^q,Q,ni,n2 3 



-Q"^q'6{u;)n,v' E 'K\9l'f\9k 

kiy , fc^ 
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Fig. 4. Diagrammatic representation of the charge density. 



Using Y.u^u\gl? [ial? + {gV'' 



2ei 



(^) ' '^'^ obtain 



(9) 



and 



V 



(3) 

Q!,j;q,Q,r2i,f!2 



(10) 



Summation over u in eq. (10) is carried out as (see Appendix), 



iDQ" 



m 



(ff] + DQ2-)^ 



and thus 



V 



(3) 



(iV(0)^)^eF iDQ"Q"'q^ 



(11) 



Using the results eqs. (7) and (11), we finally obtain the charge current at second order in J^- 
(eq. (6) ) as 



j(2)(g,Q,17i,J72) 



"Ejjfc"!'-' Q+£ ^j^iJ Q-q I ^ict'Ti W + "J^fJ + ^ilm~^ 



27] 



— '"^ 2 



(12) 



5.5 Charge density 

We calculate the charge density in the similar way. The charge density is diagrammatically 
shown in Fig. 4, and is given as 

2e/lJex^AsoW0"-imp 



p(g, Q,r2i,il2) = i 
Here, 



•^lm;q,Q,ni,n2 — / 2y,- \'^9,Q,'^,f^i,n2 ^ ^q,Q,'^,ni,n2 



'lm;q,Q,Qi,Q2 



1^=1 {fcn<^} 



l'm,3{a) I ^/»?i,3(6+c) 

:i,n2 



27r 
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^Zm,3(a) ^im,3{6+c) 
^i},Q,a;, 01,^2 + =='q,Q,cj,ni,n2 
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By using the result of Q^^Q^^^^^n^ ^"^^ O^i^^ilfS (^^- (9)), we obtain 

^ iV(0)2vr6F , 

-^«m;q,Q,f7i,Q2 ~ ^^^4 '^IV 9 > 

where A^_^ = ^^=1 E{fc„<4 7"(n:;-i5fc„_a,^)''(n:;=i5fc_„+fi,^+f2])"- The details of calcula- 
tion are presented in Appendix. After some calculation, A^ ^ is obtained as 

A°^^-l+ ^ 



Q,^ {DQ^ + in)T' 

We thus obtain ^im;q,Q,ni,Q2 + ^im;q,Q,n^,n2 = ^^g|^7''' ^i (.^'^^g2)^ - Therefore the charge 
density is 

Comparing the second term of eq. (12) (we define this term 5ja) and eq. (14), 5ja and the 
charge density p are related as 

6Uq.Q,^u^2) = iDQ''p{q,Q,VLi,^2)- (15) 

4. Result and Discussion 

Summarizing the results (eqs. (5), (12), (14) and (15)), we see that the electric current is 
written in the real space as 

3a =aiSHce«/3^ \dpS^ + [dp [s x sf + [s ^ dpSf] | - Ddo^p. (16) 

This is the main result of this paper. The first term in the curly bracket is the local contri- 
bution arising from the processes illustrated in Figs. 1 and 2. The last term daP is obtained 
by evaluating the diffusion ladder (vertex corrections) illustrated in Fig. 4. Including these 

vertex corrections is equivalent to imposing the gauge invariance by using the Ward-Takahashi 
identity.22-24 r^j^^ 

expression of the charge density (eq. (14)) is given in real space as 
p{x,t) = ^"^^"^^'^"''^ ea/j-^eQjfc j dti j dxiC{x - xi,t - ti)dkS^^^ti9jSl^ t^, 

where 

2 ^ ^i{Qt—Q-x) 

2^^I)Q2 + ij^' 
Q,n ^ 

represents the retarded diffusion propag ator satisfying {dt - DV'^)C{x,t) = 5{x)5{t). The 
consistency of our calculation is seen from the charge and current densities obtained above 
satisfying the charge conservation law, p + V • j' = 0, neglecting higher order derivatives. 

The direct (local) part of the current in eq. (16) indicates that charge current is induced if 
magnetization is dynamic, and in particular, if the magnetization has energy dissipation, SxS. 
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The current arising from the damping was found also in a Rashba spin-orbit system/^' but 
there was no contribution proportional to S in the uniform Rashba case. 

An interesting observation is that the direct contribution is expressed by the spin current 
pumped in the system. Let us define the spin current by the normal part (i.e., without spin- 
orbit interaction) as = a"' df^ip) . In the context of the spin Hall effect, proper 
definition of the spin current has been argued intensively.^^' In the present issue of the 
inverse spin Hall effect, in contrast, for the present issue of the inverse spin Hall effect, only 
the normal part of the spin current needs to be considered, since other terms related with the 
spin relaxation by the spin-orbit interaction give only higher order contribution in Ago. The 
result of the spin current is given by^^ 

jI, = c {9,5^ + ^ [d, (5 X 5)" + (5 X d,Sy] } . (17) 

We see the direct part of the charge current, eq. (16), is simply written in terms of the spin 
current, giving eq. (1). The equation for the spin current created by the spin dynamics, eq. 
(17), indicates that spin current is pumped when 1) there is either spin dynamics (S) or 
spin damping (S x S) and 2) the spin has spatial modulation (i.e., there is a finite spatial 
derivative). Applying the result to the junction system with the interface chosen as z = 0, 
perpendicular to the z axis, we may replace dzS by 5{z)S. The spin current in this case 
flows at the interface and reduces to js[z \z=o— 

cS{0){S'^ + ^^(5 x S)'^} at the interface, in 
agreement with phenomenological spin pumping theory. ^'^ 

Our weak coupling result seems very different from that for the strong coupling regime. ^^'^^ 
In ref. 28, the spin-orbit contribution to the pumped current was argued to be the scalar 
product of localized spins, ^ /5sr('^ * daS^j while in the weak coupling regime, only the 
vector product appears in the expression for the current. This situation is different from linear 
response transport properties such as the anomalous Hall effect owing to the spin structure; 
the weak and strong coupling limits in the case of the anomalous Hall effect agree when the 
spin structure slowly varies. 

The direct part of the current (eq. (16)) indicates that the term 



2t/nYT 



ceap^ I dpS^ + =^ [dp (^S X S'^ + S X dpSy j = aEf (18) 

acts as the effective field that drives the charge current, where a is the Boltzmann conductivity. 
Volovik pointed out that such a field exists in a strong coupling regime in the absence of the 
spin-orbit interaction.^*^ Volovik's electric field is written in terms of the spin variable as 
= ■ {S X daS), which is the Berry phase in space and time. This field is written 
by taking E^^^ = dtA^^^ — ^A'^^ as U(l), where A^i^^ = —0^0 cos {9 and (p are polar 
angles of the spin S and fi = 0,x,y,z). This U(l) gauge field arises by projecting the SU(2) 
gauge field of the spin onto the polarization direction. The spin gauge field in the original 
SU(2) space is given by = i(n x d^n) + lA^j^^n (n = 5/5)^°'^^ and thus Volovik's 
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Fig. 5. (Color online) Schematic illustration of the effective electric (U(l)) field arising from the 
projection of the SU(2) spin gauge field. For the strong coupling limit (adiabatic limit), the U(l) 
gauge field A^^' is that projected along the spin direction, while the dynamic inverse spin Hall 
effect is driven by the component perpendicular to the adiabatic component Aj^. 



gauge field is A^^^ = 2(n • A^). In tiie absence of the spin-orbit interaction and in the strong 
coupling limit, only this component couples to the real charge and becomes physical. This is 
understandable since only the component along the polarization direction is essential for this 
limit. In contrast, the effective electric field for the weak coupling limit with the spin-orbit 
interaction arises from the component of the SU(2) gauge field perpendicular to the adiabatic 
component. Namely, defining = 2A^ — A^^^n = n x df_in, the effective field is written as 

4^^) = ^e^^y [-dp (n X + [dpA^ + dtAj - (n x a^n)]'} . (19) 

(We note that both and n x A^ are perpendicular to n and to the adiabatic component. 
The last term is not changed since not all terms can be written using the gauge field in 
the perturbative regime). Therefore, the dynamic inverse spin Hall effect is due to a skewed 
projection of the SU(2) gauge field onto the physical U(l) field induced by the relativistic 
spin-orbit interaction. Owing to the skewed projection, the spatial coordinates and the spin 
coordinates are mixed (by the anti-symmetric tensor ea/S'y) in the effective electric field. Shibata 
and Kohno calculated the charge current induced by magnetization dynamics for the strong 
coupling limit and in the presence of the spin-orbit interaction.^^ They showed the induced 
Hall current was perpendicular to Volovik's electric field j oc Xso{n x 

Let us estimate the magnitude of the spin current and expected charge current. Assuming 
the magnetization precesses at the interface with an angle 9 and frequency /, the magnitude 
of the spin current density is evaluated from eq. (17) as 

b«l = l + 6^sm^ / A/m . 



37r2 V ^ y V 

For a disordered metal with ~ 0.1 coupled weakly to the magnetization (^ = 0.01), the 
spin current density is estimated to be |js| = 1 x 10^'^ K/rv? for a frequency of / = 10 GHz^ 
and /cp = 1 X 10^^ 1/m^ if ^ 0(0. 1). We estimate the spin-charge conversion efficiency as 
oiSH = 1 X (10^'^ ~ 10~2) using the dimensionless spin-orbit coupling fcpAso = 0.01 0.1 
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for various materials. We therefore expect the magnitude of the output charge current to be 



5. Conclusions 

In conclusion, we presented gauge-invariant theory for the dynamic inverse spin Hall effect 
satisfying charge conservation. Assuming metallic systems, the case of symmetric spin-orbit 
interaction due to random impurities was considered. We showed that charge current is induced 
directly by the effective electric field, which is proportional to the local spin current, but not 
by the diffusive process. The effective field is explained as a skewed U(l) projection owing to 
the relativistic spin-orbit interaction. 
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Appendix: Derivation of the diffusion factor 

We present, in this section, details of the calculation of the infinite series summation due 
to the vertex corrections.^^' In the calculation, we see 



|j| ~ 1 X (10^ ~ 10^) A/m^. 




oo 




Above equations are rewritten as 



A° - =A" - il A° - (A° - )^ + f A° - )^ H I = A" - II -h A' 

- =A° - +(A° ~)'^ + (A° -\ = A° - 11 A° - 1 




(A-1) 



where we define 
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Expanding A" - and - with respect to Q, f2 in slowly varying limit, we obtain as 

Since we obtain A^ ^ = ^— ^ from eq. (A-1), Aq and Ag are given as 

~„ iPQ- 



A°Q = -1 + 



1 



{DQ^ + in)T 
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